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Abstract 

We consider a class of parametrically forced Hamiltonian systems with one-and- 
a-half degrees of freedom and study the stability of the dynamics when the frequency 
of the forcing is relatively high or low. We show that, provided the frequency of 
the forcing is sufficiently high, KAM theorem may be applied even when the forcing 
amplitude is far away from the perturbation regime. A similar result is obtained 
for sufficiently low frequency forcing, but in that case we need the amplitude of the 
forcing to be not too large; however we are still able to consider amplitudes of the 
forcing which are outside of the perturbation regime. Our results are illustrated by 
means of numerical simulations for the system of a forced cubic oscillator. In addition, 
we find numerically that the dynamics are stable even when the forcing amplitude 
is very large (beyond the range of validity of the analytical results), provided the 
frequency of the forcing is taken correspondingly low. 

Keywords: Nonlinear oseillators, high driving frequeney, low driving frequeney, 
perturbation theory, KAM theorem averaging. 
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1 Introduction 

Forced systems can display unexpected behaviour in the extreme cases in which the os¬ 
cillation of the forcing is either very slow or very fast. It is well known that, in such 
situations, the dynamics can be strikingly different with respect to the case in which the 
forcing period is comparable with that of the unperturbed motions mEH[121E13ED]. A 
classical example is provided by the pendulum with oscillating support ESElIiaElElls]. 
The standard technique used to attack the problem is the averaging method; however, 
a rigorous implementation is quite not trivial, as it involves dealing with small divisor 
problems. Indeed, the main idea underlying the averaging is that the forced system, in 
suitable coordinates, can be considered as a perturbation of an integrable system, so that 
KAM-like arguments apply [ 26 l [I] . 


I 


For a periodically forced system, if the amplitude of the forcing is small, then the 
KAM theorem yields that most of the unperturbed tori persist. On the basis of heuristic 
arguments, one expects that something of the same kind still occurs even when the ampli¬ 
tude of the forcing is not small, provided its frequency is large or small enough. Such an 
expectation is ultimately based on averaging arguments. However, to reduce the problem 
to results available in the literature requires some analysis. Since this has not been done 
explicitly, as far as we are aware of, we think it may be of interest to discuss how the 
argument works. 

The case of large frequency can be reduced quite easily to Neishdatd’s averaging theo¬ 
rem. However, at a certain point, we shall need an additional close-to-identity transforma¬ 
tion which strongly relies one the one-dimensionality of the system (more precisely on its 
integrability). The case of low frequency is related to the problem of boundedness of the 
solutions of forced systems in phase space jniEniiiziEiiES]. However, to prove bounded¬ 
ness one has to prove the existence of confining KAM tori far away from the origin. Hence 
only the asymptotic behaviour of the potential really counts, and in general one needs a 
condition on the growth of the potential at infinity - besides smoothness conditions, see 
[231 EH- On the contrary, to study the existence of KAM tori in a fixed region of phase 
space of a forced system, one needs information about the potential in that region: this 
explains why the assumptions we shall require on the potential are stronger, as they are 
not just asymptotic properties. 

For clarity purposes, in this paper we focus on a case study, the forced cubic oscillator, 
which has been extensively investigated in the literature all along as a paradigmatic model 
ISlIlllsllig. We shall prove that most of the phase space is filled by KAM invariant tori 
when the period of the forcing is sufficiently large or sufficiently small. In particular we find 
that this happens for values of the forcing amplitude far beyond the perturbation regime. 
Furthermore, we provide numerical results which strongly suggest that the amplitude can 
be allowed to be still larger than the analytical bounds. We leave as an open problem how 
to improve the analytical bounds, so as to give full explanation of the numerical findings. 
We briefly discuss at the end how to relax the assumptions on the potential and extend 
the results to more general systems. 

A cubic oscillator subject to a periodic driving force is described by the Hamiltonian 

2 4 

H{y,x,t) = ^ + {1 + (1.1) 

where (^, x) G R^, t denotes the time, the driving force / is a 27r-periodic analytic function 
of its argument, with zero average and \\f\\oo = 1; /x G R and cj G R are, respectively, 
the amplitude and the frequency of the driving force. The phase space for the system is 
R2 X T^, with = R/(27r/6c;)Z and the corresponding equations of motion are 

'x = y, 

< y =-{1 + (1.2) 

4 = r 

where the dot denotes derivative with respect to t. 
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For /i = 0 the system is integrable, so that the full phase space is filled by invariant 
tori, with t increasing linearly in time with frequency 1 and (x, y) moving on a closed orbit 
C in the plane with frequency Q. depending on the initial data. 

Fix arbitrarily ^^2 > > 0 and consider, for the unperturbed cubic oscillator, the 

closed, concentric orbits Ci and C 2 which run with frequencies and respectively. 
Call Vq C the bounded region enclosed between Ci and C 2 , and set V — VqX T^. The 
reasons behind the choice of the domain V are the following: we want to fix a bounded 
region in phase space so as to estimate the relative measure of the persisting tori (the 
curve C 2 is a natural boundary for such a region) and at any rate one is forced to exclude 
a small region around the origin, where chaotic dynamics may become dominant (the role 
of the internal curve Ci is just to cut off such a small region since the beginning). 

For fi small enough, we can apply KAM theorem [1], so as to obtain the persistence of 
most of the invariant tori inside P, independently of the value of the frequency u (and of 
the value of as well: can be taken arbitrarily small). 

Theorem 1. For y small enough the set of persisting invariant tori in V for the system 
with Hamiltonian 0 leaves out a set with relative measure 

In this note we want to show that the similar results still hold when removing the 
condition of smallness on /i, provided the frequency uu is either large enough or small 
enough. This will be discussed, respectively, in Section and Section In particular, by 
increasing cj from 0 to infinity, a double transition regularity-chaos-regularity is expected to 
occur. We support such an expectation by providing in Section [^numerical results, which 
also give evidence that the regularity regime extends to wider ranges of the parameters for 
which the analytical results do not apply. Finally in Section we discuss how to extend 
our analysis to more general systems. 

2 High frequency regime 

We consider first the case of uu large and we setcj — l/s^ with 6: small. We can formulate our 
result for such a case in a more general context. In act ion-angle variables, the Hamiltonian 
0 becomes 

H— Hq (A) “h iJjH\ (a, tr, cjt), (2.1) 

where 

= J cn\Ta)f{t), (2.2) 

with cna := cn((a, 1/v^) and T := 4iF(l/v^)/27r, where cii{a,k) and K{k) denote re¬ 
spectively the cosine-amplitude function and the complete elliptic integral of the first kind 
with elliptic modulus k. 

The corresponding Hamilton equations are 

(a = noiA) + iidAHi{A,a,Lot), 

M = -ij,daHi{A,a,Lot), 
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where rio(^) OaHo^A) = . In terms of the variables (A^aA) the domain V 

defined in Section [^becomes 2) = 2to x T x where T = Ti and 2to = {A G R+ : Ai < 
1^1 < ^ 2 }, with Ai = nlT^/3 and GI 2 = n^T^/S. 


More generally we shall consider Hamiltonians of the form (2.1) and make the following 
assumptions on Hq and Hi — trivially satisfied in the case (2.2). 


Hypothesis 1. Assume the Hamiltonian funetion (2.1) to he real-analytie in a domain 
0 :=2lo X T X where 2lo is an open subset o/R. 


Hypothesis 2. Assume A 1 -^ i2o(^) dAHQ{A) to he a loeal diffeomorphism on SIq. 

f At 

Hypothesis 3. Assume {Hi{A^ cr, •)) := / — Hi{A^ cr, t) — 0^ ^{A^ cr) G 2to x T. 

Jt 2vr 


In (2.1) rescale time t ^ r — out. Then the equations of motion (2.3) become 

{ a = £:i 2 o(^) + £iadAHi{A,a,r), 

A = -£iadaHi{A,a,T), 


(2.4) 


where now the dot denotes derivative with respect to time r. 

We can apply Neishtadt’s averaging theorem [SUES! to cast the system into the form 


a — £ (rio(A) + /adAVei^^ + l^dARe^A^ <^5 ' 7 ”)) 5 
A = -£ {jidaVe^A, a) + ladaReiA ^)) , 


where and are suitable analytic functions, with VAA^a) = a, •)) + O(^) 

and Rs an exponentially small remainder, that is \Rs\ < (7exp(—c/s) for some positive 
constants c and C. The change of coordinates is canonical and s-close to the identity. In 
order not to overwhelm the notations, we denote the new variables with the same letters 
as the old ones. 

By Hypothesis]^ the average of Hi vanishes, hence I 4 is a correction of order s to i^o- 
So one can perform a further close-to-identity change of coordinates which leads to the 
equations 

(a = ene{A) + siJ.dARe{A,a,T), 

\a^ -eiJ,daRs{A,a,T), 

where i2£(A) = i2o(^) + 0(s) and R is still exponentially small (and again we still denote 
by (A^a) the transformed coordinates). The corresponding Hamiltonian is 

H{A, a,£)=£ (HAA) + yiRe{A, a, r)) . (2.7) 

The overall change of coordinates leading to ( |2.7| ) is close to the identity within 0{£) 
and hence, up to a region with measure 0(6:), the domain 0 is transformed into a region 
enclosed between two KAM invariant tori. By studying the Hamiltonian H{A^a^£)/£ we 
see that we can apply once more KAM theorem and conclude that most of the unperturbed 
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tori for the Hamiltonian Hs{A) persist when the perturbation Rs{A^a^r) is switched on. 
Since now the perturbation is exponentially small, the relative measure of the tori which 
are destroyed is exponentially small in To go back to the original coordinates, we have 
to scale back time. So we obtain the following result. 


Theorem 2. Consider the system with Hamiltonian (2.1) and assume Hypotheses^to 
For any value of ji, for uo large enough the domain 2) is filled by KAM invariant tori, 
up to a region whose relative measure is 0{1 /uj). Apart from a thin region elose to the 
boundary, the invariant tori leave out a region with measure exponentially small in l/uo. 


Note that Hypothesis is needed here, contrary to Theorem to ensure that the 
averaged system is integrable (such a condition is automatically satisfied for /i small, 
without any assumption on Hi). 


3 Low frequency regime 


Now consider o> with (jj = e. We can reason as in m (see also da HU Emm). Fix 
|/x| < 1, so that 1 + > 0. We rewrite the equations of motion 0 as 


(3.1) 


\x = y, 

\y = -a{£t)x^, a{t) ■-1 + 

Then the argument proceeds through the following steps. 

First, through a time-dependent canonical change of coordinates {x,y) (A, a), with 

/q4\2/3 

j (a(£t))“^/®cn(ra), y = -^—j (a(£t))^/®sn (Ta) dn (To), (3.2) 

where sn{a) and dn((a) are the sine-amplitude and delta-amplitude functions with modulus 
k — respectively, one writes (3.1) as the Hamilton equations corresponding to the 

Hamiltonian 


T-L(A, a, t) = T-Lo(A, a, t) + cn (Ta) sn (Ta) dn (Ta) -^ 7 ^, 

6 T a[et) 

no{A,a,t) 


- {a{et)y 


(3.3a) 

(3.3b) 


where b{t) = d{t) and T := 4iF(l/v^)/27r. Next, in order to eliminate the dependence on 
time in the leading term one makes the change of coordinates {A, a,t) 1 -^ (p, q, s), with 


q^t, 


s = a, 


p = l-L{A,a,t), 
which leads to the Hamiltonian 

^(PW,s) = AQ{p,£q)+£Ai{p,£q,s), 

T {4:pfA 


(3.4) 


Ao{p,q) = 


3(a(Q))V4’ 
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(3.5a) 

(3.5b) 






for some function Ai of order 1 in e. Then one may perform a further change of variables 
(p, g, s) (J, (/p, s) into action-angle variables for Ao^ so yielding the Hamiltonian 


B{J, if, s) = Bq{J) + eBi{J,e(p, s), (3.6a) 

Bo{J) = {njf\ (3.6b) 

for a suitable function >81 of order 1 in 6 :. 

Finally we integrate the Hamiltonian equations corresponding to ( |3.6| ) between s = 0 
and s = 27 t. Denote by the solution; then, defining = £(p{s) and setting 

= (J(27r),'0(27r)) and (J,'0) = (J( 0 ),' 0 ( 0 )), we obtain the twist map 


J' = J + £2f(J,V’), 

+ ef2o( J) + e'G(J, ^P), f2o( J) = 


for suitable analytic functions F and G. 

Therefore we can apply Moser’s twist theorem and conclude that any invariant curve 
with Diophantine rotation number persists for £: small enough. For fixed 6 :, the relative 
measure of the persisting curves in a given region of the cylinder is O(v^) [HI EO]. Note 
that, in order to prove just boundedness of the solutions, it would be enough to prove 
the existence of an invariant circle of constant type (as sometimes done in the literature, 
see for instance HH; see also the comments in m)- On the contrary, to prove that the 
persisting tori have large measure for 6 : small, a milder Diophantine condition is required; 
one could even allow a Bryuno condition on the rotation number [ 8 ], as done in [ii> but 
this would not increase appreciably the measure of the invariant curves. 


Coming back to the original coordinates we obtain the existence of a large measure set 
of invariant tori for the continuous flow which has the twist map (3.7) as Poincare section 
at times multiples of 27r. We can summarise the discussion by the following statement. 


Theorem 3. Consider the system with Hamiltonian O) and fix fi ^ ( — 1,1). For uj 
small enough the domain 2) is filled by KAM invariant tori, up to a region whose relative 
measure is 


4 Numerical results 


Numerically we may illustrate the scenarios considered in the previous sections as well 
as provide some insight into cases which are not covered by the analysis. We study the 
system (1.1) with f{t) — cost. In particular we consider three situations, |/i| < 1, |/i| > 1 
and |/i| = 1 , with both high and low frequency forcing. 

To check the stability of the dynamics we take 10 000 pseudo-random initial conditions 
within a square [—2,2] x [—2,2] from the phase plane {x,y). The chosen numerical inte¬ 
gration method is a Stormer-Verlet scheme with variable step size. The Stormer-Verlet 
method is a second order symplectic scheme, details of which may be found in m- 
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After an initial transient period, the trajectories are checked to ascertain how their 
asymptotic behaviour has changed with respect to the trajectories of the unperturbed 
system with the same initial conditions. If most of the orbits have remained close to the 
corresponding orbits of the unperturbed system, then we say that the system is “stable”. 
This is expected to occur when the system is well within the KAM regime: the majority 
of the unperturbed tori persist slightly deformed, so that every orbit either lies on a torus 
or is trapped between two surviving tori. However it is possible that the trajectories of 
the perturbed system do not remain close, but are still bounded. This can happen as 
we are moving out of the KAM regime: most of the tori are destroyed, with a few of 
them still existing and undergoing much larger deformations. We refer to such a case by 
saying that the system is “bounded”. Numerically it is difficult to classify trajectories 
as unbounded, as a trajectory which appears unbounded may be bounded within a very 
large region. Therefore, pragmatically, we class the trajectories as unbounded once their 
amplitude exceeds 30 in either the x oi y direction, and class the system as “unbounded” 
if any trajectory is found to be so. When this happens, nearly all (if not all) KAM tori are 
expected to be destroyed, at least in the region investigated, otherwise any of them would 
confine the orbits inside. We note, however, that, even though KAM theory no longer 
applies in this case, one can still have invariant curves of a different kind, such as cantori 
(see for instance [28]), so that it is quite possible for some trajectories to be unbounded 
whilst others remain bounded (and possibly not too far from those of the unperturbed 
system). 

First we consider the scenario where |/i| < 1, in particular we choose /i = 0.8, 0.9 
and 0.95; some numerical results are shown in Table For uu sufficiently large or small, 
as proved in the previous sections, the dynamics are stable, whilst between these two 
extremes the system loses stability. 


UJ 

II 

o 

bo 

y = 0.9 

y = 0.95 

0.05 

Stable 

Stable 

Stable 

0.1 

Stable 

Stable 

Bounded 

0.2 

Stable 

Bounded 

Bounded 

0.3 

Bounded 

Bounded 

Unbounded 

0.4 

Bounded 

Unbounded 

Unbounded 

0.5 

Bounded 

Unbounded 

Unbounded 

0.7 

Bounded 

Unbounded 

Unbounded 

0.8 

Unbounded 

Unbounded 

Unbounded 

12.0 

Unbounded 

Unbounded 

Unbounded 

13.0 

Stable 

Unbounded 

Unbounded 

14.0 

Stable 

Stable 

Stable 


Table 1: Stability and boundedness results for choices of |/i| < 1 and various values of uj. We 
class the system as “stable” if most of the orbits remain close to those of the unperturbed system, 
“bounded” if none of the orbits exceed 30 in either direction and “unbounded” if it is not bounded, 
i.e. at least one orbit is unbounded; see text for details. 
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In Figure we show some example orbits corresponding to the initial conditions 
{x^y) = (1,1) with /i = 0 in Figure a) and y = 0.8 in (b), (c) and (d). It may be 
seen that the perturbed orbits remain close to the unperturbed system for suitable uj. 
With fi getting close to 1, one has to take oj increasingly small or increasingly large for 
the system to be stable. 




(a) (b) 




(c) 


(d) 


Figure 1: Orbits for the system (1.2) with initial conditions {x,y) = (1,1). In Figure (a) y = 0. 
In Figures (b), (c) and (d) y = 0.8 and uj = 0.0001, 0.2 and 14, respectively. 




(a) (b) 


Figure 2: Poincare maps for the system (1.1) with y 
the break up of the KAM tori. The frequency is u; = 


= 0.8 and various initial conditions, showing 
0.2 and 0.4 in (a) and (b), respectively. 































The results in Table [T show that for ji — 0.8 the system becomes unstable when 
uj > 0.2. We see in Figure ^b) and (c) that as u increases from 0.2 to 0.4, many of the 
tori are broken and cantori appear, separating the few persisting KAM tori and creating 
larger and larger gaps. However, as quite a few tori still persist in the region investigated, 
the orbits remain bounded. Increasing cj further causes more of the KAM tori to break 
up and for u = 0.8 the orbits are no longer bounded in the region [—30,30] x [—30,30]. 
Analogous considerations apply to the other values of fi: the closer the amplitude fi is to 
1 , the more extreme u must be for the system to be stable. 


UJ 

/i = 1 

/i = 1.2 

(1 = 2 

/i = 5 

0.0001 

Stable 

Stable 

Stable 

Stable 

0.0002 

Stable 

Stable 

Bounded 

Unbounded 

0.0003 

Bounded 

Unbounded 

Unbounded 

Unbounded 

0.0004 

Bounded 

Unbounded 

Unbounded 

Unbounded 

0.0005 

Bounded 

Unbounded 

Unbounded 

Unbounded 

0.0010 

Unbounded 

Unbounded 

Unbounded 

Unbounded 

13 

Unbounded 

Unbounded 

Unbounded 

Unbounded 

14 

Stable 

Stable 

Unbounded 

Unbounded 

16 

Stable 

Stable 

Unbounded 

Unbounded 

18 

Stable 

Stable 

Stable 

Unbounded 

20 

Stable 

Stable 

Stable 

Stable 


Table 2: Stability and boundedness results for choices of |/i| > 1 and various values of uj. The 
system is classed as either “stable” or “bounded” or “unbounded” as explained in the caption of 
Figure 


For |/i| > 1 the analysis in the previous sections can only be applied when the system 
undergoes high frequency forcing. Numerically we find that the system is also stable with 
low frequency forcing, however u must be taken considerably smaller than the cases where 
|/i| < 1. This is not true for high frequency forcing, where similar orders of lj (compared 
with the cases where |/i| < 1) are sufficient for the dynamics to become stable. Some 
numerical results are presented in Table In Figure we show some example orbits with 

/i > 1. 

Similarly to the case where |/i| < 1 it is evident that the perturbed orbits remain close 
to that of the unperturbed system provided u is suitably chosen. In Figure [^f) we see 
that, although the system is classed as unbounded for uj = lA and /x = 5, it is still possible 
to find initial conditions for which the orbit remains close to the unperturbed system. 


5 Conclusions 

More generally the Hamiltonian of a forced cubic oscillator is 

H{y,x,t) = + {a + (5.1) 
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However the writing (1.1) is not restrictive, since we can reduce (5.1) to that form by 
rescaling both variables x and y and redefining the parameter y. 






Figure 3: Orbits for the system (0 with initial conditions (x,^) = (1,1). In Figures (a),(b) and 
(c) uj = 0.0001 and in Figures (d), (e) and (f) uj = 14. The amplitude is /i = 1, 2 and 5 in Figures 
(a, d), (b, e) and (c, f), respectively. 


Other generalisations can be easily envisaged. For instance any potential V(x) yielding 
closed orbits in a region encircling the origin can be considered. In particular one can take 
a potential instead of x^/4: the unperturbed system is still integrable, so that 

the analysis of Section applies immediately. Also the discussion in Section can be 
easily adapted to cover such a case; we refer to nanzi for details. Also, less regularity 
can be required for the driving force. Finally, one could consider quasi-periodically forced 
systems, as in [22], in the case in which all components of the frequency vectors are small 
or large. 

Coming back to our system (O), the condition that / has zero average could be 
relaxed. In fact, for Theoremj^to hold, what we really need is that l + /i(/(-)) > 0, so that 
the averaged system is integrable. On the other hand. Theorem requires 1 + yf(’) > 0. 
Thus, if (/(•)) = 0 and ||/||oo = C this excludes the case /i = 1. However, if the average of 
1 + fifiut) is positive, one can argue that the potential remains positive for most of time, 
so one can conjecture that the condition 1 + y{f{-)) > 0 might be sufficient in the low 
frequency regime as well. As discussed in Section we have found numerically that also 
in such a case the orbits are stable if the forcing frequency is sufficiently low. It would be 
interesting to investigate the issue in more detail by means of analysis. 
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